CREATING WAVE-FOCUSING MATERIALS 
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ABSTRACT 

Basic ideas for creating wave-focusing 
materials by injecting small particles in a 
given material are described. 

The number of small particles to 
be injected around any point is calcu- 
lated. Inverse scattering problem with fixed 
wavenumber and fixed incident direction of 
the plane acoustic wave is formulated and 
solved. 

1. INTRODUCTION 

This paper contains the results presented 
at the author's plenary talk at the IPDO- 
2007 symposium on inverse problems, design 
and optimization. A method for creating ma- 
terials with a desired refraction coefficient is 
given. This coefficient, in particular, may be 
chosen so that the new material has a desired 
wave-focusing property (see also [1]-[12]). 

We want to investigate the following 
problem. Let D be a bounded domain filled 
with a material whose properties are known, 
for example, a homogeneous material with 
known speed propagation of sound waves. 
Can one inject into D small particles in such 
a way that the resulting new material would 
have some desirable wave-focusing proper- 
ties? For example, is it possible to create in 
this way a material that scatters an incident 
plane wave in a desired solid angle? 

There is a large engineering and physi- 
cal literature on creating "smart" materials. 
Photonic crystals, quantum dots, coating, 
are some key words. However, there seems to 
be no prior work which deals with the ques- 
tion posed above. 

In this paper this question is studied 
rigorously: the "smallness" of the particles 
is specified, the number of these particles 
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around the point x is specified, the role of 
the shapes of these particles is explained, and 
the notion of wave-focusing is made precise. 
The basic results of this paper are: 
1) It is proved that the injection of a suit- 
able number of small, acoustically soft parti- 
cles in a given bounded region, filled by some 
material with known properties, allows one 
to create a new material such that its scatter- 
ing amplitude is arbitrarily close to a given 
scattering amplitude. 

2) A method is given to calculate the den- 
sity N(x) of small particles, to be injected 
in a unit volume around a point x 6 D, in 
order that the new material has the scatter- 
ing amplitude close to the desired scattering 
amplitude. 

3) For the first time the problem of finding 
a compactly supported potential q{x) which 
generates the scattering amplitude A(f3) := 
A(f3,a,k), approximating an arbitrary fixed 
given function /(/?) <E L 2 (S 2 ) with any 
desired accuracy is formulated and solved. 
Here the wavenumber k > and the inci- 
dent direction a € S 2 are fixed, S 2 is the 
unit sphere in R 3 . 

In Section 2 detailed statements of the 
problems are given. In Section 3 the inverse 
scattering problem with fixed wavenumber 
k > and fixed incident direction a G S 2 
is discussed. 

In Section 4 the ill-posedness of the above 
problems is discussed. 

For engineers the paper gives a "recipe" 
for creating materials with a desired refrac- 
tion coefficient. No such receipes were given 
earlier, to the author's knowledge, although 
there were many papers (see book [13] and 
references therein) in which an effect of em- 
bedding small spheres or ellipsoids into a ho- 
mogeneous material on the effective dielec- 
tric and magnetic properties of the material 
was discussed, and the new material was a 
homogenized material. Homogenization in 



the literature was studied mostly for elliptic 
positive-definite operators, while we study 
wave propagation and the corresponding op- 
erator is not positive-definite. 

Another principal difference between this 
paper and the earlier published results con- 
sists in the statement of the problem: we are 
not considering the distribution of small par- 
ticles, embedded into the material, as uni- 
form or random, but solve a design problem 
of creating materials with a desired refrac- 
tion coefficient. 

Moreover, we give (see also [10]-[12]) a 
precise recipe for creating such material for 
practically arbitrary refraction coefficient. 

The two technological (engineering) prob- 
lems to be solved for this recipe to be prac- 
tically implemented, can be formulated pre- 
cisely as well: 

a) How does one practically embed in 
a given material many small particles given 
the number of the particles per unit volume 
around every point of the original material? 

b) How does one prepare practically 
small particles with a desired boundary 
impedance? 

In this paper we consider small parti- 
cles with the Dirichlet boundary condition, 
corresponding to acoustically soft particles. 
Problem b) will arise if one considers the 
small particles on the boundary of which an 
impedance boundary condition is imposed. 
Varying the boundary impedance as a func- 
tion of positions of particles, one can create 
refraction coefficients with the desired ab- 
sorption properties. This was discussed in 
more detail in [12]. 

2. STATEMENT OF THE PROBLEM 
AND SOME RESULTS 

Consider a bounded domain DcK 3 with 
a smooth boundary S. The scattering of a 
plane wave on this domain is described by 
the equations 

[V 2 + k 2 nl(x)]u = in D, k = const > 0, 

(1) 

u = e lka ' x + Ao(0, a,k) — +o (-) , 

r W (2) 

r := \x\ — > oo, (3 := — . 

r 

Here a £ S 2 is given, A ((3, a, k) is the scat- 
tering amplitude, k > is fixed through- 
out the paper, n^x) > is a given func- 
tion, the refraction coefficient, tIq(x) = 1 in 
D' := R 3 \D, Uq(x) is piecewise-continuous. 
The function n^(x) describes the material 
properties of the region D. Problem (l)-(2) 
has a unique solution u £ Hf oc {R 3 ), where 



Hf oc (R 3 ) is the Sobolev space. Equation (1) 
can be written as a Schrodinger equation 

[V 2 +k 2 -q Q (x)]u = 0, q (x) := k 2 [l-n 2 (x)}, 

qo = in D'. Suppose that M small 
acoustically soft particles (bodies) D m are 
injected into domain D. Smallness means 
that ka << 1, where a is the characteris- 
tic size of the small bodies. One may define 
a := | maxi< m <M diamD m . Assume that 
the boundaries S m of D m are uniformly Lip- 
schitz, i.e., the Lipschitz constant docs not 
depend on m. Acoustically soft means that 
u\s m — 0, where u can be interpreted as 
acoustic pressure. 

The scattering problem in the new region 
can be formulated as follows: 

M 

[V 2 + k 2 n 2 {x)]U = in R 3 \ (J D m , (3) 

m— 1 

U\ Sm =0, 1 < m < M, (4) 

e ikr / ]\ 

U = u+Am(P, a, k) ho-, 

r W (5) 

x a 
r 

where u solves (l)-(2). One can replace (5) 
by the following: 

ikr / -i \ 

U = e lka ' x + A(P, a,k) + o - 

r W (6) 

x a 
r — > oo, — = p. 

r 

From (2), (5) and (6) one gets A = Aq + Am- 

Since n^x) is known, A is known. 

Problem (3)-(5) has a unique solution. 

Problem 1 is to show that one can dis- 
tribute sufficiently large number M of small 
particles in D in such a way that A((3, a) := 
A(P, a, k), k > is fixed, will approximate in 
L 2 (S 2 ) an arbitrary fixed scattering ampli- 
tude /(/?, a) with any desired accuracy. We 
consider also the following problem. 
Problem 2: Can one distribute small parti- 
cles in D so that the resulting new material 
would have the scattering amplitude A{[3) 
(k > and a £ S 2 are both now fixed), 
which approximates in L 2 (S 2 ) an arbitrary 
given function /(/?) £ L 2 (S 2 ) with any de- 
sired accuracy? 

If A(/3, a) := A(/3, a, k) {k > is fixed) 
is a scattering amplitude, known for all 
(3, a £ 5 2 , then one can find the unique, 
corresponding to A{(3,a), potential q(x) by 
the Ramm's method [1], [2]. This method 
gives a stable approximation of q even in 



the case when noisy data Ag(f3, a) are given, 
su P/3, Q es2 \A(0, a) - As(/3, a\ < S. 

If q(x) is found from A(/3, a), then we 
define p(x) := q — qo(x) and prove that 
p(x) = N(x)Co, where N(x) is the number 
of small particles per unit volume around a 
point x, i. e. , the spatial density of the num- 
ber of the particles, and Co is the electrical 
capacitance of a small conductor of the same 
shape as the particle. Here we assume that 
all the small particles are identical, but this 
assumption can be dropped (see [3], [4] ). 

Let us summarize: If one injects small 
particles with the spatial density N(x) = 
where p(x) := q(x) —q a (x), then the re- 
sulting new material will have practically the 
desired scattering amplitude A(/3, a), corre- 
sponding to a potential q(x). 

This gives a solution to Problem 1. 
Note that the scattering amplitude A(/3, a), 
V/3, a G 5 2 , corresponding to a real- valued 
potential q G L 2 (D) 1 determines q uniquely. 
Since qo(x) is known, the function 



N(x) = 



q(x) - q Q (x) 
Co 



gives an approximate solution to Problem 
1. This solution is not exact because we did 
not pass to the limit 
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but took just sufficiently small identical par- 
ticles with Co being electrical capacitance of 
the perfect conductor of the same shape as a 
single particle. 

The number of small particles per unit 
volume is 0(-^). Their volume per unit vol- 
ume of the original medium is 0(§r)- This 
quantity tends to zero as | — > 0. The ca- 
pacitance per unit volume is O(Jf). Thus, 
the limit of the ratio Jj is finite and non- 
zero, while the limit of the relative volume 
of the injected small particles is zero because 
^3 tends to zero. 

3. INVERSE SCATTERING WITH 
FIXED k AND a 

Our solution to Problem 2 is based on 
the idea used in Section 2 in solving Prob- 
lem 1. Given /(/3), we find q(x) G L 2 (D), 
such that the corresponding to q scattering 
amplitude A q (f3) (k > and a G S 2 are 
fixed) approximates with a desired ac- 
curacy: - A q ((3)\\ L 2 {D2) < e, where 
e > is an a priori given small number. If 
such a q is found, then N(x) = g ^ c ^°^ as 
in Section 2. The principally novel problem 



is finding q from /(/3) and e. This problem 
has many solutions, as we prove. However, a 
priori it is not at all clear if this problem has 
a solution. Let us outline our solution to this 
problem. First, recall the well-known exact 
formula for the scattering amplitude: 



AM 



J e ikf) ' x q{x)u(x)dx 1 



where u(x) is the scattering solution, 
[V 2 + k 2 - q{x)]u = in i? 3 , 

U = Uq + A q ((3)—— + o(-), T ► OO, 

P = f, u := e ika ' x , a G S 2 and k > 
are fixed, u{x) = u(x,a,k). Denote h{x) = 
q{x)u{x). Then 



-ik/3-x 



' x h(x)dx, 



where A q is the scattering amplitude, corre- 
sponding to the potential q. Given f((3) and 
e > 0, however small, one can find (many) h 
such that 

H/0 9 ) + ^ / e-^h(x)dx\\ L2(s2) < e. 

(7) 

The function h can be found, for example, 
as a linear combination h n = Y^j=i c i¥'j( x )i 
whercj^} is a basis of L 2 (D). If n is suffi- 
ciently large and Cj are found from the min- 
imization problem 



n ,. 



-ik(3-x n 



<Pj{x)dx\[i l 2i S 2\ = min, 



then (7) holds. In [4] another, analytical, 
solution to (7) is given. 

If h = h £ {x) is found, then q{x) can be 
found from the nonlinear equation h = qu. 
This equation for q is nonlinear because the 
scattering solution u = u(x;q) depends non- 
linearly on q. One has 



u(x) 



u (x)- / g(x,y)q(y)u(y)dy, (8) 
Jd 



where u = e lka ' x and g := 47r|x _ y| 
h(x) 



q(x) :-- 



u (x) - J D g(x,y)h(y)dy' 



(9) 



If the right side of (9) is an L 2 (D) function, 
then (9) solves our inverse scattering prob- 
lem. Indeed, define 

u(x) := u (x) - / g(x,y)h(y)dy. 
Jd 



This u solves (8) with q defined in (9). The 
scattering amplitude 



ik/3-x 



h(x)dx. 



By (7) one has 

\\f((3) - A q ((3)\\Li( S i) <e. 

So, Problem 2 is solved if (9) defines an 
L 2 (D) function. This, for example, is the 
case if 



inf \uo(x) 



D 



g( x ,y) h (y)dy\ > °- 



If formula (9) defines a non-integrable func- 
tion due to possible zero sets of the function 



ip(x) 



:= mo - / g(x,y)h(y)dy, 
Jd 



then, as we prove, a suitable small pertur- 
bation hs of h in i 2 (_D)-norm will lead to a 
function 

• =- | e L'W. 



where 



1p5 



= u - g(x,y)h s (y)dy, 
Jd 



and \\h — hsW^iD) < 8. Since the function 
— J D e~ lk @' x hs(x)dx differs a little from 
the function — j-J D e~' lk ^' x h(x)dx, condi- 
tion (7) is satisfied if h is replaced by hg and 
e by, for example, 2e. Therefore 



Ng(x) :-- 



qgjx) - qojx) 
Co 



solves Problem 2 in the case when (9) is not 
an I/ 2 (£>)-function. 

Let us explain how to choose hg. First, 
without loss of generality one can assume h 
to be analytic in D = D [J S, because an- 
alytic functions (even polynomials if D is 
bounded) are dense in L 2 (D). If h is ana- 
lytic, so is ip(x) in D. Therefore the null set 
of ip, 

M := {x : ip(x) = 0, x G £>}, 
is generically a line defined by two equations 
ipi := Retp = 0, ip2 '■= Imtp = 0. 

Let 

Ms ■= {x : \ip\ < 5, x G D}, 

and D s := D\Mg. Generically |V^| > c > 
on M and, by continuity, this inequality 
holds in Ms (possibly with a different c > 0) . 



Small perturbation of h leads to these generic 
assumptions. 
Define 



hs 



h in Ds, 
in Mg, 

<7 A ■— J U <>-J D 9(x,y)h s (y)dy d ' 

| in Ms. 



(10) 



Then q s G L 2 (D) and 



N s (x) r- 



qsjx) - qojx) 

C 



solves Problem 2. 

Let us check that qs G L 2 (D). We prove 
more: qs G L°°(D). It is sufficient to check 
that 

inf \ips(x)\ = 

xGDs 

inf \u (x) - / g(x,y)h(y)dy\ > cS > 0, 
x ^ D s Jd s 

(11) 

because qs = in Ms by the definition. 

Choose the origin on M and make a 
change of variables 

si = ipi(x), s 2 = tp2(x), s 3 = x 3 . (12) 

The Jacobian of this transformation of vari- 
ables is non-singular because V-01 and VV>2 
are linearly independent on TV and in Mg. We 
have 

maxdJI + IJ" 1 !) <c, 

xeM s 

c > stands for a generic constant indepen- 
dent of 5. Let us check that 



We have 



where 



\^s{x)\ > cS in D s . 



\M*)\>\*K*)\-i(8), 



1(6) = f \ 9 {x,y)h{y)\dy. 

If x G Ds, then > 8, and 

\Mx)\> 5-1(5). 

Moreover, max^u |ft| < M and, using the 
new variables (12), one gets 



1(5) <c[ dpp [ 
Jp<S Jo 



ds 3 



because the region Ms can be described by 
the inequalities 

p 2 = s\ + sj < S 2 , 0<s 3 <l, 



and we have used the estimate |J 1 | < c in 
Ms- Integral 

1(5) = 0(6 2 \\n5\) as 5^0. 

Thus 

IV^I >6-0(S 2 \\nS\) > cS 

with some constant c 6 (0,1). This justi- 
fies our method for solving Problem 2 in 
the case when formula (9) does not yield 
q e L 2 (D). 

4. ILL-POSEDNESS OF PROBLEMS 

1 AND 2 

Both Problems 1 and 2 arc ill-posed. 
Since the ill-poscdness of Problem 1 has 

been discussed in great detail in [1], [2], we 
discuss only the ill-posedness of Problem 2. 

In Problem 2 one has to find h, given 
/ and e > 0, so that (7) holds. This is an 
ill-posed problem, similar to solving the first 
kind integral equation 

Bh := J e~ ikl3x h{x)dx = /(/?) 

for h. If this equation is solvable for a given 
/, it may be not solvable if / is replaced by a 
slightly perturbed function f$. If e > in (7) 
is small and / is not in the range of B, then 
||/i||i,2(£)) is large. This leads to large maxi- 
mal values of the corresponding q. Therefore 
any numerical method for solving Problem 

2 has to use a regularization procedure. In 
[5] one can find some numerical results re- 
lated to Problem 2 and a description of the 
regularization procedure which was used. 
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